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Traditional spintronics relies on spin transport by charge carriers, such as electrons in semiconduc- 
tor crystals. This brings several complications: the Pauli principle prevents the carriers from moving 
with the same speed; Coulomb repulsion leads to rapid dephasing of electron flows. Spin-optronics 
is a valuable alternative to traditional spintronics. In spin-optronic devices the spin currents are 
carried by electrically neutral bosonic quasi-particles: excitons or exciton-polaritons. They can form 
highly coherent quantum liquids and carry spins over macroscopic distances. The price to pay is a 
finite life-time of the bosonic spin carriers. We present the theory of exciton ballistic spin transport 
which may be applied to a range of systems where bosonic spin transport has been reported, in 
particular, to indirect excitons in coupled GaAs/AlGaAs quantum wells. We describe the effect 
of spin-orbit interaction of electrons and holes on the exciton spin, account for the Zeeman effect 
induced by external magnetic fields, long range and short range exchange splittings of the exciton 
resonances. We also consider exciton transport in the non-linear regime and discuss the definitions 
of exciton spin current, polarization current and spin conductivity. 

PACS numbers: 71.35.-y, 03.75.Kk, 03.75.Mn, 73.63.Hs, 78.55.Cr 



I. INTRODUCTION 



Excitons are electrically neutral and have finite life- 
times. These are two obstacles which prevented the de- 
velopment of excitonic spintronics, or spin-optronics un- 
til the beginning of this century. How possibly one can 
speak about the current, which is nearly independent on 
the applied bias, and whose amplitude changes dramat- 
ically with distance and time? - was a typical question 
by the sceptics of exciton-based spintronics. While elec- 
trons and holes have been considered as perfectly valid 
spin carriers, and exotic effects like the spin Hall effect 1 
have enthusiastically been searched for— ~ , the spin cur- 
rents carried by excitons 7 -"— and exciton-polariton o 10 i 1 1 
over tens or even hundreds of micrometers received only 
limited attention by the scientific community. There ex- 
isted a huge imbalance of theoretical works on fermionic 
and bosonic spin transport. This is changing now. More 
and more researchers work on bosonic spin current prob- 
lems 1 ^—, and this research is paid for by beautiful recent 
discoverie o 1 1 1 18 ~ — . To summarize tens of publications in 
one sentence: bosonic systems bring new quantum coher- 
ent effect s 21 ' 22 to the physics of spin transport. Stimula- 
tion^— and amplification 2 ^ of spin currents are possi- 
ble in exciton and exciton-polariton Bose gases. Besides 
bosonic effects, exciton spin transport has another impor- 
tant specific feature: it is dissipative by its nature, as the 
spin carriers have a finite (and short for excitons in reg- 
ular materials) life-time. In continuous wave optical ex- 
periments stationary spin textures can appear: excitons 
are injected in the structure, they propagate ballistically 
or diffusively, and eventually disappear by radiative re- 
combination. Their polarization properties and spin are 
inherited by photons emitted by the system, which is 



why the polarization patterns observed in near field pho- 
toluminescence experiments directly characterize exciton 
spin currents in the plane of the structure. 

The goal of this work is to define what the exciton 
spin, magnetization and polarization currents are, and 
to explain how they can be described within the most 
frequently used spin density matrix (DM) approach and 
mean-fiel d 27 ' 28 approximation. We consider a specific 
system, namely a planar zinc-blend semiconductor struc- 
ture containing quantum wells, where excitons can be 
formed. This choice is motivated by recent experimen- 
tal discoveries in GaAs/AlGaAs based coupled quantum 
wells. We limit the scope of this paper to heavy-hole 
excitons, however, our approach can be easily extended 
to light-hole excitons or excitons in quantum wells of a 
different symmetry. We do not speak here about the 
large variety of recent experimental results and applica- 
tion of the formalism presented here to the description of 
one particular experiment, as this would make this pa- 
per too long and too specific. For a direct comparison 
of theoretical simulations with the experimental data we 
address the reader to Ref. H|. We are confident that the 
approaches formulated here are suitable for the descrip- 
tion of a great variety of excitonic spin effects in quantum 
wells. 

In Section II we introduce the spin density matrix 
DM formalism accounting for the different mechanisms 
of spin re-orientation and the relation to electron and 
hole spin currents. In Section III we present numerical 
results obtained within the spin density matrix formalism 
and analyze them. In Section HVl we study the non-linear 
spin dynamics of propagating excitons using the Gross- 
Pitaevskii (GP) equations. The next three sections of 
the paper are devoted devoted to exciton spin currents 
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and polarization currents, 
are given in Section VIII. 



Conclusions are perspectives 



II. THE SPIN MATRIX FORMALISM FOR 
PROPAGATING EXCITONS 

In zinc-blend semiconductor quantum wells (e.g. in the 
most popular GaAs/AlGaAs system), the lowest energy 
exciton states are formed by electrons with spin projec- 
tions on the structure axis of +1/2 and -1/2 and heavy 
holes whose quasi-spin (sum of spin and orbital momen- 
tum) projection to the structure axis is +3/2 or -3/2. 
Consequently, the exciton spin defined as the sum of the 
electron spin and heavy-hole quasi-spin may have one 
of four projections on the structure axis: -2,-1, +1, +222. 
These states are usually nearly degenerate, while there 
may be some splitting between them due to the short 
and long-range exchange interactions. 

It is important to note that the present formalism ad- 
dresses the spin part of the exciton wavefunction, which 
is a product of electron and hole spin functions. For ex- 
ample, the probability to find the exciton in the spin state 
+1 is given by a product of probabilities to find an elec- 
tron in the spin state -1/2 and the heavy hole in the spin 
state +3/2. The four component exciton wave-function 



* = (* +1 ,*_i,# +2! #- 2 ) 
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where \& e + i and ^e.-J are the components of the elec- 
tron spinor wavefunction; ^> h + 3 and ^ h _s are the com- 
ponents of the heavy-hole spinor wavefunction. 

To describe the dynamics of the system we will first 
define the Hamiltonian (Section III Ap describing the dif- 
ferent physical mechanisms of spin evolution. We then 
introduce the spin density matrix DM (Section III Bp for 
the description of exciton spin states. We relate its com- 
ponents to the observable Stokes' vectors of light emitted 
by bright excitons and use the Liouvillc equation to de- 
scribe its evolution in time. From the density matrix 
DM one can describe the polarization state of excitons 
and thus exciton spin currents. It is also instructive to 
consider the relationship with electron and hole spin cur- 
rents (Section III Cp . 



Hamiltonian 



discussed in the Section ITVI and neglecting magnetic field 
effect on center-of-mass motion and internal structure of 
excitor*22r— . We consider excitons propagating ballisti- 
cally in the plane of a quantum well. We shall character- 
ize them by a fixed wave-vector, k ex . We represent the 
full exciton Hamiltonian as a sum of three parts de- 
scribing the spin-orbit and Zeeman effects on electrons, 
T-L% x , holes, THLg X , and the exchange induced splittings of 
exciton states. 
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Dresselhaus and Zeeman terms 



In order to build the 4x4 matrix Hamiltonian for exci- 
tons, we start with simpler 2x2 Hamiltonians describing 
the spin-orbit and Zeeman effects for electrons and holes. 

The electron Hamiltonian in the basis of (+1/2,-1/2) 
spin states is: 

1 



He = Pe (k e ,x<Jx ~ k e . y Uy) - -p e /i B B5-. 



(3) 



Here g e is the electron g-factor, \xb is the Bohr magne- 
ton, B is a magnetic field, a is the Pauli matrix vector, 
and p e is the Dresselhaus constant describing spin-orbit 
interactions of electrons. The Pauli matrix operators are: 
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(4) 

Rewriting Eq. ([3]) and retaining only z-component of the 
magnetic field, which corresponds to the Faraday geom- 
etry, one can obtain: 



-He = 



-\g<,\x B B p e (k, 

Pe {k e ,x - ik e ,y) \g e ^B B 

-\ge\i B B p e k e e^ 
Pete-* \ge\x B B 



(5) 



where <j) is the angle between k ox and the chosen x-axis. 
The exciton Hamiltonian needs to be written in the basis 
of (+1,-1, +2,-2) exciton spin states, which correspond to 
(-1/2, +1/2, +1/2, -1/2) electron spin states. The elec- 
tron spin-flip couples +1 and +2 states as well as -1 and 
-2 states. For each of these two couples of states we apply 
the Hamiltonian [Eq. ([5])]. which results in the following 
electronic contribution to the 4x4 exciton Hamiltonian: 
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Here we derive the exciton Hamiltonian in the basis 
of -2,-1, +1, +2 states, accounting for the spin-orbit in- 
teraction (Dresselhaus and Rashba effects)^, long- and 
short-range exchange interactions^ and Zeeman effect, 
but neglecting exciton-exciton interactions, which will be 



2. Rashba terms 

We note that another possible spin-orbit contribution 
to the Hamiltonian may come from the Rashba effect, 
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which takes place in biased quantum wells. The Rashba 
term to be added in Eq. is a e (a x k e y — u y k e x ), where 
a e is a constant proportional to the Rashba field. The 
contribution of the Rashba term to the electron Hamil- 
tonian, in the basis of (+1/2,-1/2) electron spin states, 
can be re-writtcn: 



HL — 





-ia e k e e l 



ia e k e e 




(7) 



Using the same procedure as for Dresselhaus terms, this 
gives an additional contribution to the exciton Hamilto- 
nian, in the basis of (-1/2, +1/2, +1/2, -1/2) electron 
spin states: 
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(8) 

Unless stated explicitly, wc will for simplicity omit the 
Rashba terms in the rest of this paper and consider only 
the Dresselhaus terms. Note that there are no linear in 
wave-vector Rashba terms for heavy holes in zinc-blend 
quantum wells grown along the (OOl)-axis. 



3. Heavy-hole contribution (Faraday geometry) 

The heavy hole contribution to the Hamiltonian can 
be calculated with similar reasoning to the electron 
case. The heavy-hole Hamiltonian written in the basis 
of (+3/2,-3/2) states is: 

Hh = Ph (kh,x&x + kh,yd- y ) - -9h^BBa z . (9) 

Here gh is the heavy-hole g-factor and f3h is the Dressel- 
haus constant for heavy holes^. 
obtain: 



Re- writing Eq. (|9]), we 
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(10) 



The hole spin-flip couples +1 and -2 states as well as -1 
and +2 states. For each of these two couples of states we 
apply the Hamiltonian (fTOj) , which results in the following 
hole contribution to the 4x4 exciton Hamiltonian: 
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4- In-plane magnetic field (Voight geometry) 



quantum wells. Suppose that the field is applied in the 
x-direction. In the electron and hole basis the Zeeman 
Hamiltonian is in this case: 



He 
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-^9e,h^BBa x 







— 29e,hHBB 




(12) 

Note that the hole g-factor in plane of the quantum well 
is different from the g-factor in Faraday configuration, 
in general. This maps into the (+1, —1, +2, —2) exciton 
basis as a Zeeman Hamiltonian of the form: 



Hz = - 



(j-bB 



g e g h 

g h g e 

9e gh 

gh g e 



(13) 



5. Exchange terms 

Besides the contributions from electron and hole spin 
orbit interactions and Zeeman splitting, there may be a 
purely excitonic contribution to the Hamiltonian, which 
is composed from the Hamiltonian for bright excitons 
written in the basis (+1,-1): 



Hb = E b I - 5b<J x = 



E b -S b 
S b E b 



(14) 



and the Hamiltonian for dark excitons written in the ba- 
sis (+2,-2): 



Hd = EdI - Sd<T x 



E d S d 
S d E d 



(15) 



where / is the identity matrix. The terms with 8 b and Sd 
describe the splittings of bright and dark states polarized 
along x and y axes in the plane of the structure due to the 
long-range exchange interaction. E b — Ed is the splitting 
between bright (+1 and -1) and dark (+2 and -2) exciton 
states due to the short-range exchange interaction. In 
microcavities, this splitting is additionally enhanced due 
to the vacuum field Rabi splitting of exciton-polariton 
modes formed by bright excitons and a confined optical 
mode of the cavity^. 

The origin of Eqs. (fT4")) and (|T5j) can be easily seen from 
the exciton Hamiltonian written in the basis of linear x 
and y polarizations. For example, for the bright excitons: 



Hxy 



E b -S b 
E b +, 

H b = C 1 HxyC, 



(16) 
(17) 



where: 



If the magnetic field is applied in the plane, it splits 
electron and hole states polarized in the plane of the 
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are the transformation matrices from the linear to circu- 
lar polarization basis and vice versa^I. The same reason- 
ing can be applied to the dark excitons as well. 

The sum of Hamiltonians lib and Tid, written in the 
4x4 exciton spin basis is: 



Now, the full exciton Hamiltonian can be written as: 



nje 
ft, 
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For the translational motion of an exciton as a whole 
particle the exciton momentum is given by ~P ex = 
(m e + nth) ~v e x , where m e and are in-plane effective 
masses of an electron and of a heavy hole, respectively; 
v ex is the exciton velocity. Having in mind that the ex- 
citon translational momentum is a sum of electron and 
hole translational momenta given by P c .h = mejiV e .ft, 
where v e and v/j are the electron and hole velocity, re- 
spectively, one can easily see that v/, = v e = v ex . 
Having in mind that P ex = Hk ex and P e h = ^-ehi 



we have k e 



k h + k e , with k e 



and 



k ex . Thus, H ex depends on the exciton cen- 
ter of mass wave-vector k ea; and on the angle <fi between 
this angle and one of the structure axes (e.g. (lOO)-axis). 

It should be noted that in this consideration the wave- 
vectors k ea; , k/j, k e are related to the translational motion 
of the exciton as a whole particle, hole and electron as 
its constituents. The wave-vector of relative motion of 
the electron and hole "inside" the exciton is zero on av- 
erage but may be important for each given moment of 
time. Recently, the effect of relative electron-hole motion 
on the spin-orbit effects of excitons has been analyzed by 
Vishnevsky et . Their analysis confirms the presence 
of linear in k ea; spin orbit terms in the exciton Hamilto- 
nian introduced above. 



B. Spin Density Matrix 

Having constructed the Hamiltonian for excitons prop- 
agating with a wavevector k ox , we now consider the de- 
scription of their spin state. We shall use the spin density 
matrix, p = |*) (*|, where * = (vf( +1 , vfi_ 1 , v]/ +2j vfi_ 2 ) 
are the components of the exciton wavefunction projected 
onto the four spin states, (|*+i) , , |*+2) , |*-2))- 



1. Relation to Stokes' vectors and polarization degrees of 
light 

The elements of the upper left quarter of the density 
matrix DM are linked to the intensity of light emitted 
by bright exciton states, /, and to the components of the 
Stokes' vector, S x , S y and S z of the emitted light: 



I 



Pll 



S z 



P12 — S x — iS y , 

P21 = S x + iSy, 

I 



P12 



s z 



(21) 

(22) 
(23) 

(24) 



These expressions can be summarized more succinctly 
using the Pauli matrices as: 



Pu Pl2 
Pl\ P22 



I 



-I + S.a, 



(25) 



where S = (S x , S y , S z ) is the Stokes' vector and we recall 
that I is the identity matrix. 

Often when studying the polarization structure of 
fields with non-uniform intensity, it is useful to compare 
the polarization degrees of emitted light, which can be 
given by normalizing the Stokes' vectors to the light in- 
tensity. The circular polarization degree is: 



2S Z 



Pc = 



Pll - P22 



I Pll + P22 

The horizontal- vertical linear polarization degree is: 



(26) 



Pi 



2S X 



Pl2 + P21 
Pll + P22 



(27) 



The linear polarization degree measured in the diagonal 
axes (also referred to as a diagonal polarization degree) 
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is given by: 



25 
~1 



_ ^Pl2 — P21 
Pll + P22 



(28) 



2. Liouville equation 



The dynamics of the density matrix DM is given by 
the quantum Liouville equation: 



ih 



dp 



njtot ' 



(29) 



where the Hamiltonian is composed from the electron, 
hole and exciton contributions given by Eqs. ©, (fTTj) 
and (|T9| (considering the Faraday magnetic field config- 
uration). 

So far, we have neglected all relaxation or scattering 
processes in the system. The commonly used way to ac- 
count for these processes is through the introduction of 
a phcnomenological Lindblad supcroperator to the Liou- 
ville equation: 



ih = 



dp 
It 



njtot ~ 
rl *ex ' r 



L(p), 



(30) 



where the Lindblad superoperator is introduced as: 



L (p) = ih 



pu/n pu/n pu/tc pu/t c 

P2l/U P22/U P23/t c P2i/T c 

P3l/l~c P32/t c P33/Td P3i/Td 

Pai/t c . PviItc Pa3/t c i Paa/t c i 



(31) 



Tb is the bright exciton decoherence time, Td is the dark 
exciton decoherence time and r c is the characteristic de- 
coherence time of processes between dark and bright ex- 
citons. Note that dissipation is crucial in the description 
of exciton spin currents. In particular, within this for- 
malism, the current conservation and flux conservation 
should be completed in the continuity equation by exci- 
ton generation and decay. 



The formalism described so far, in sections III Al and 
this section pi Bj) . has been successfully applied in the 
description of spin transport in gases of cold excitons in 
coupled GaAs/AlGaAs quantum wells^. In this work, 
cold excitons are generated within localized spots and 
then fly away ballistically in radial directions. The ele- 
ments of the density matrix, pij, are dependent on the 
distance from the excitation spot r = v ex t at time t 
and the polar angle, <fi. The propagation speed v ex = 
h\i ex / (m e + rrih). In this work, when solving the Liou- 
ville equation introduced above, we shall refer to the ex- 
perimental configuration of Ref. H|. In particular, this 
implies a specific choice of the initial conditions for Eq. 
(|30|) : we shall assume that at zero time excitons are 
not moving. We shall assume that they populate the 
cigenstates of the exciton Hamiltonian taken with 
h ex = following a thermal distribution with a tempera- 
ture T. We shall neglect all relaxation processes assum- 
ing L (p) = 0. 

C. Electron and hole spin currents 

We shall normalize exciton, electron and heavy-hole 
functions to unity, namely: 



** , 1 * P 4.1 + 1 * P 1 



■ e,+ _ 



3 + 



1 (32) 
1 (33) 

1 (34) 



Now, the exciton spin density matrix DM is given by: 



l*X*l 
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*!j_ 2 *_i *l 2 *_i 
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(35) 



which is represented in terms of electron and hole wave- 
functions as: 
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136] 



This representation allows us to obtain useful links be- 
tween the elements of exciton, electron and hole density 



matrices, in particular: 



Pe = |*e) (* e | 



Ph = 



$ 1 

e 



P22 + P33 P24 + P31 
Pl3 + PA2 Pll + PAA 



(37) 



ft,+- 



ft, + 5 



ft, 



Pll + P33 P14 + P32 
P23 + PA^ P22 + PAA 



(38) 
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We know that the components of electron and hole den- 
sity matrices are linked with the projections of electron 
and hole spins as: 



P< 
Ph 



S R 



s, 



e,x * 
1 



IS, 



Sh,x + iSh,y 2 — &h,z 



(39) 
(40) 



where, for the heavy hole we have assigned spin +1/2 to 
the state +3/2 and spin -1/2 to the state -3/2 accounting 
for the orbital momentum of these states of +1 and -1, 
respectively. 

The z-component of the spin polarization carried by 
electrons can now be expressed as: 



Se, 



(P22 + P33 - Pll ~ Pii) /2. 



(41) 



Similarly, the z-component of the spin polarization car- 
ried by holes can now be expressed as: 



Sh,z = (Pll + P33 ~ P22 - Pte) /2. 



(42) 



The in-plane component of electron and hole spins can 
be extracted from the off-diagonal elements of the den- 
sity matrix. Namely, the x-component of electron spin is 
given by: 

S e ,x = (Pl3+/>31+/>24 + P42)/2, (43) 

while the x-component of the hole spin is given by: 

S h ,x = (pi4 + P23 + P32 + P4i) /2. (44) 
The y-component of electron spin is given by: 

S e ,y = * {~Pl3 + P31 + P2i - Pi2) A (45) 

while the y-component of the hole spin is given by: 

Sh,y = i {Pli - P23 + P32 - P4l) /2, (46) 
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FIG. 1: (color online) (a) Dispersion of the excitonic states 
calculated using the set of parameters from Table U (green 
lines) and reduced set of parameters with Et — Ed — 0, St = 
0, Ph = (red symbols). Black arrows indicate the values 
of the wavevectors used in (c)-(e) k ex = 15.3^m _1 and (f) 
k ex = 5/im -1 . Double-ended arrows indicate the energies of 
oscillations between the eigenstates which appear in spatial 
polarization patterns, (b) Boltzmann distribution at T=0.1K 
chosen as initial condition at the source for density matrix, (c) 
Linear polarization degree along x-axis calculated with the 
simplified set of parameters, (d) Same as (c) but for the full 
set of parameters from the Table [T] (e) Circular polarization 
degree with parameters from Table [I] (f) Same but at k ex = 
5/im" 1 . The source area was taken circular with a radius of 
4/im. 



TABLE I: Parameters for numerical calculations. 



Electron mass 


m e 


0.07 


mo 


Heavy hole mass 


m h 


0.16 


mo 


Electron Dresselhaus coupling 




2.7 




Heavy-hole Dresselhaus coupling 


Ph 


0.92 


^ieV/im 


Bright exciton XY splitting 


5 b 


0.5 


^eV 


Dark exciton XY splitting 


6 d 


-13 


AieV 


Bright-dark exciton splitting E 


— E d 


5 





III. NUMERICAL RESULTS IN THE DENSITY 
MATRIX FORMALISM 

Figure Q] shows the numerical results obtained within 
the density matrix formalism for a model system with 
the same parameters as those of coupled double quantum 
wells studied in Ref. [H The parameters arc summarized 
in Table U 

The dispersion of bright and dark exciton modes is 
shown in Fig. [TJa) by green solid lines. The momentum 
has been chosen along the x direction, but the anisotropy 
of the band structure remains small. The numerically cal- 
culated linear and circular polarization degrees are shown 
in Figs. QJd) an d ( e )> respectively. In FigQJd), four lobes 
are unambiguously observable in the pattern of the linear 
and diagonal (not shown in figures) polarizations. This 
pattern is a consequence of the Dresselhaus spin-orbit 



coupling for electrons, and it is characteristic of the cho- 
sen initial state: the thermal distribution of excitons over 
four split eigen-states having zero in-plane wave-vector. 
The exciton distribution function chosen in this calcu- 
lation as initial condition at the source is shown in Fig. 
QJb). It corresponds to the temperature of 0.1 K. At such 
a low temperature, the initial state is dominated by the 
lowest lying state, which is a linearly polarized dark ex- 
citon. This choice of initial conditions leads to a variety 
of polarization patterns observed in experiment^. 

In order to reveal the mechanism of formation of the 
polarization patterns it is instructive to consider a sim- 
plified version of the Hamiltonian [20] In the absence 
of magnetic field, zero Dresselhaus effect for holes, zero 
splitting of bright excitons and zero splitting Eh — Ed be- 
tween bright and dark excitons this Hamiltonian can be 
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rewritten as: 



H = -S d 



£e~ # o 
^ 
1 
£e"^ 1 



(47) 



where £ = —(3 e k e /8d- This Hamiltonian has an ana- 
lytical solution, whose eigenvalues are given by: E = 
±i(5d ± 5<5d\A + 4£ 2 , corresponding to two branches at 
low energy oc ±£ 2 and two branches at high energy 
oc ±(1 + £ 2 ) for £ <C 1. The corresponding bands are 
shown by red squares in Fig. [TJa) and the eigenvectors, 
starting from the lowest energy and taking 5d < 0, can 
be approximated for small £ by: 
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As mentioned already, at low temperatures fegT <C |6<i|, 
the lowest energy state with no momentum is given by 
[0, 0, —1, 1] and is a linearly polarized dark exciton. After 
an acceleration due to dipolc repulsion, this initial state 
is no longer an eigenstate and oscillates as a function of 
time between the two eigenstates of the Hamiltonian with 
which the initial state is not orthogonal, namely the first 
and the third eigenstates listed in Eqs. [48] Among these 
two states, only the third gives a significant contribution 
to the observed polarization as it is essentially "bright" 
(has large projections to +1 and —1 exciton states). The 
linear and diagonal polarizations originating from this 
state are readily given by — cos(2</>) and — sin(20) re- 
spectively, which reproduces the essential features of the 
numerical results for the linear polarization pattern, as 
one can see comparing the images calculated with the re- 
duced Hamiltonian (Fig. [TJc)) and the full Hamiltonian 
(Fig. (Hd)). In the particular case considered here (ini- 
tial state formed essentially by linearly polarized dark 
excitons) the Dresselhaus spin-orbit term for electrons 
leads to formation of the linear polarization vortex: the 
polarization plane is always perpendicular to the wave 
vector direction. The linear polarization vortex has been 
observed experimentally by High et al^ 

We should emphasize that with a proper reordering of 
the basis vectors, the simplified Hamiltonian 27] is analyt- 
ically equivalent to the Hamiltonian of bilaycr graphenc. 
In its simplest expression, the Hamiltonian of bilayer 
graphene can be written as^S: 
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(49) 



where vf is the velocity, k the amplitude of the mo- 
mentum, 4> the angle between the momentum and the 



x axis, and t± is the main coupling term between the 
two graphene layers, which is given by hopping between 
two carbon atoms that are superimposed. The four coef- 
ficients of the associated wavevectors correspond to the 
probability amplitudes on the two independent sublat- 
tices of the two graphene layers. By permutation of the 
basis vectors:l — s- 1, 2 -> 3, 3 — > 4, 4 — s- 2, the Hamilto- 
nian is rewritten as: 
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(50) 



which is equivalent to Eq. 37J with t± = —5d and 
v F k/t± = £. 

The 4x4 Hamiltonian of bilayer graphene is often 
restricted to the subspace of the two low energy bands. 
From Eq. 1481 the corresponding eigenstates are: 



(51) 



We can then define a pseudospin vector, which represents 
the relative phase between the two components of the 
wavevectors. From Eq. 1511 the pseudospin rotates two 
times when the particle wavevector undergoes one full 
rotation. 

Adapting the terminology used for graphene, the four 
lobes observed in Fig. [ljc,d) would be nothing else than 
the fingerprint of the "pseudospin" rotation in the ex- 
citon system. In the context of the exciton system un- 
der consideration, the phase <j> corresponds to the angle 
between the exciton polarization vector and the chosen 
x-axis of the structure. 

The build up of circular polarization requires intro- 
duction in the model of the splitting between linearly 
polarized bright exciton states 5b- This splitting acts as 
an effective magnetic field applied to the Stokes vector 
of light emitted by bright excitons, S = (S x , S v , S z ). In 
particular, if the x-polarized exciton state has a lower 
energy than the y-polarizcd exciton state, it creates an 
effective magnetic field in the ^-direction which rotates 
the Stokes vector in the yz plane. This converts the diag- 
onal polarization to the circular polarization and leads to 
the appearance of right- and left-circularly polarized sec- 
tors in the polarization map of exciton emission. Separa- 
tion of spins due to lincar-to-circular polarization conver- 
sion is known in exciton-polariton systems as the Optical 
spin Hall effect. This effect was theoretically predicted 
in Ref. [ID and experimentally observed in polaritonioi^ 
and excitonic systems^. Note that the circular polariza- 
tion pattern is strongly sensitive to the chosen exciton 
wave- vector, which governs the energies of the four in- 
volved eigen-states. Figure QJe) was obtained with the 
wave vector k ex = 15.3/im -1 corresponding to the cross- 
ing point of the dispersion branches associated to the 
first and second exciton eigen-states (see Fig[TJa)). To 



compare. Fig. QJf) shows the circular polarization pat- 
tern calculated with k ex = 5/im _1 . One can see that the 
four lobes pattern of circular polarization is washed out, 
because the lowest energy state of the system remains 
essentially dark. Rapid oscillations that show up, have 
the period determined by the splitting between the two 
lowest states. 

Finally, let us underline that the dispersion curves 
shown in Fig. HJa) do not take into account the kinetic 
energy of excitons. The kinetic energy would shift all 

B 2 fc 2 

curves up by K = „ . Note that this does not affect 



electron 



hole 



the splittings between exciton eigenstates and would not 
affect the spin dynamics of excitons. In the spin density 
formalism developed above we assign to all excitons the 
same kinetic energy, K. In a realistic system, the kinetic 
energy may be spread, in which case averaging of the 
obtained polarization patterns over k ex may be needed. 
This averaging would smoothen the fast oscillations seen 
on the images (c,d,e,f). In the next section devoted to the 
non-linear spin dynamics we will fully take into account 
the kinetic energy of propagating excitons. 

The distribution of in-plane projections of electron and 
hole spins for the same choice of parameters as above is 
shown in Figs.[^a-h). The left panels show electron spins 
and the right panels show the hole spins. The directions 
of spins are shown by arrows, while the length of each ar- 
row is proportional to the value of the spin polarization. 
The upper panels, Figs. [2Ia,b), show the spin distribu- 
tions in the absence of a magnetic field. In this case 
the electron and hole spins are oriented along the effec- 
tive Dresselhaus fields which are oriented differently for 
electrons and heavy holes, as we have discussed in the 
previous section. The magnetic field strongly changes 
the spin distribution in real space. The spin textures be- 
come strongly anisotropic in the case of in-plane (x- or y- 
oricntcd) magnetic field. Note, that the in-plane isotropy 
in the system is broken by the splitting between x- and 
y- polarized exciton states, which is why switching of 
the magnetic field between x and y axes strongly affects 
the distribution of electron and hole spins. It should be 
noted also that electron and hole in-plane spin textures 
can hardly be observed directly in optical experiments. 
However, they can be deduced from fitting the exciton 
polarization maps, e.g., using the formalism described 
above. 

Figure |3] shows how the magnetic field affects spatial 
patterns of linear (a-d) and circular (e-h) polarization. 
Switching the magnetic field orientation between the x-, 
y- and z-axes one can dramatically affect the polarization 
patterns. Having in mind that the exciton polarization 
patterns can be directly observed in near-field photolu- 
minescence experiments, fitting of these patterns to the 
experimental data would allow extracting the Dressel- 
haus constants and exciton exchange splittings, which, 
in turn, allow to restore electron and hole spin textures^. 

Figure BJa-d) illustrates a peculiar regime where the 
Dresselhaus fields for electrons and holes are taken to be 
zero and there is no magnetic field applied, but electrons 
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FIG. 2: (color online) Spatial distribution of electron and hole 
spin components in the plane of the quantum well structure, 
calculated at zero magnetic field and at B — 5 T at three 
different orientations. The parameters are given in Table [T] 
at k ea; = 15.3/im -1 .The source area was taken circular with 
a radius of 4fj,m. 



are subjected to the Rashba field (the Rashba field for 
heavy holes is zero). Figure BJa) shows the electron spin 
distribution in space, where the spins are clearly aligned 
along the Rashba field force lines. Interestingly, the hole 
spins become polarized as well, while no effective field 
acts on them (Figure Bib)). This is an illustration of the 
exciton effect: bound in excitons by Coulomb interaction 
and subject to the exchange induced exciton effects, the 
holes acquire in-plane spin polarization. The non-zero 
spin polarization of heavy holes is possible due to the ex- 
citon exchange effects. Figure EJc,d) shows the resulting 
linear and circular exciton polarization patterns. One can 
see that the Rashba effect induces polarization patterns 
strongly different from those induced by the Dresselhaus 
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FIG. 3: (color online) Spatial distribution of exciton linear 
polarization ((a)-(d), length of the bars maps the intensity in 
arbitrary units) and color maps of the exciton circular polar- 
ization degree ((e)-(f)), same color code as in Fig. [TJ. Mag- 
netic field B = in (a), (e), and B — 5T along z in (b),(f), 
x in (c),(g), y in (d),(h). Parameters are given in Tabled] 
k £ z = 15.3/an -1 . The source area was taken circular with a 
radius of 4/xm. 

effect, which is why from the shape of polarization pat- 
terns one can conclude on the nature of spin-orbit cou- 
pling in the system. 



IV. NON-LINEAR SPIN DYNAMICS OF 
PROPAGATING EXCITONS AND 
EXCITON-POLARITONS 

In the previous section we operated with a spin den- 
sity matrix which is very convenient for the description of 
partially coherent and partially polarized exciton gases. 
The quantum Liouville equation (|30[) is a very efficient 
tool for the description of effects linear in the exciton 
density. On the other hand, one cannot straightforwardly 
incorporate non-linear interaction terms in this equation. 
The treatment of non-linear effects in a partially coherent 
system is a non-trivial task. Much simpler is the treat- 
ment of non-linear effects in a perfectly coherent system, 




-20 20 -20 20 



x (|jm) 

FIG. 4: (color online) Spatial distribution of electron (a) and 
hole (b) in-plane spin component calculated at B=0 in the 
absence of the Dresselhaus field /3 e = 0; f3h — 0, but including 
Rashba field a e = 2.7^eV/^m;. Other parameters are given in 
Tabled] k ea; = 15.3/im" 1 . The source area was taken circular 
with a radius of 4fim. Corresponding patterns of linear (c) 
and circular (d) exciton polarization degree are shown with 
the same color code as in Fig. [T] 

such as a condensate at zero temperature. In this case, 
the ensemble of excitons can be described by a single 
4-component wave function = $f-i, ^+2, ^-2) 

. The linear dynamics of this wave-function for ballis- 
tically propagating excitons having a wavevector ]z ex is 
described by the Schrodinger equation: 

ihj t \y)=U\*), (52) 

where the Hamiltonian is the same as in Eq. (f2"U)) . This 
equation represents a set of four coupled linear differen- 
tial equations for four exciton spin components. Non- 
linear effects lead to the condensate evolution in real and 
reciprocal space. 

From now on we shall consider the exciton spin dynam- 
ics in real space (2D), so that the wave function ^ will be- 
come coordinate-dependent and will not be restricted to 
one single value of k ex . The non- linear interaction terms 
for multi-component exciton gases are introduced and 
discussed in detail in Ref. 0. Here we expand Eq. (|52|) 
by introducing the kinetic energy (to describe the real 
space dynamics) and the interaction terms. On the other 
hand, we neglect the magnetic field, for simplicity. This 
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results in a system of four non-linear Schrodinger or GP 
equation a 40 " 42 : 

fr 2 V 2 (3 e m e ( . n 
ih—— = * + i H fcz - f +2 

, Phrrih 



k x - iky) *_ 2 + ai |*+i| *+i 



a 2 |*-i| 2 *+i + a 3 |*+2| 2 *+i 
a 4 |^_2| 2 *+i + W*I 1 * +2 *_ 2 , (53) 
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h 2 v 2 



/3 e m e 
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a 2 |*+i| 2 *-i + « 3 |*-2| 2 *-l 

|* +2 | 2 *_1 + W^* +1 ^+2^-2 



k x +ik y ) * +2 +ai |*_i| *_i 



(54) 



ih- 



dt 



h 2 v 2 

2m ex 
Phm h 



(3 e m e 



+2 



(k x + iky) *+i 



+ ^ (fc - 4) *-i + ai |* +2 | 2 * +2 

+ a 2 |*_ 2 | 2 * +2 + a 3 |*+i| 2 * +2 

+ a 4 |*_i| 2 * +2 + W r * , l 2 * + i1'_i, (55) 



ih 



ft 2 V 2 


2772 err 



fc^ - iky ) 



a 2 |* 



+2| 



a 3 |*. 



+ a 4 |*+ir*-2 + W r *+2*+i*-i- (56) 

Here k x ^ y = —tV x ,yi m ex = tn e + f^ft- To make this 
system more compact we have omitted the terms de- 
scribing exchange induced exciton splittings given by the 
Hamiltonian (f^Uj) . We do not discuss here the nature 
and value of the interaction constants ai : 2,3,4 and W. In 
the system of indirect excitons in coupled GaAs/AlGaAs 
quantum wells, as a zeroth approximation, one can take 
cii = q 2 = «3 = a 4 . Note also that in microcavities, 
where the lower exciton-polariton mode is strongly de- 
coupled from dark excitons, the dark exciton states may 
be almost empty at low temperatures. If this is the case, 
the spin dynamics of the exciton-polariton condensate 
is given by the first two of the four GP equations (|53l 
and I54p with 0^4 = W = 0. The remaining constants 
ai j2 have been widely discussed in literature. We address 
the interested reader to the experimental and theoretical 
study of these constants by Vladimirova et al 4 ^. 

The GP equations are widely used for the description 
of coherent propagation of exciton-polaritons in micro- 
cavities^ 4 -. They allow for the studying of interesting 



topology effects such as: quantum vortices 4 ^— ; half- 
quantum vortices^—; bright 5 ^ and dar k 49 ' 55 " — solitons. 

The polarization of light emitted by an exciton or 
exciton-polariton condensate can be obtained as: 
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(57) 



(58) 



(59) 



These expressions easily follow from the definition of the 
spin density matrix. 

A significant limitation of the GP equations as a theo- 
retical tool is that they assume a coherent state of the sys- 
tem. If one is interested in the spin structure of the zero 
temperature ground state of excitons in a Bose-Einstein 
condensate 4 ^, then this assumption is fulfilled by defini- 
tion. However, in real systems there is an incomplete co- 
herence that, strictly speaking, requires a description of 
statistical mixtures, perhaps involving density matrices. 
Furthermore, Eqs. (|53M56p have been written assuming 
an infinite lifetime for the particles (be they excitons or 
exciton-polaritons) , which is never the case of real non- 
equilibrium systems. Often pumping and radiative decay 
terms are introduced into Eqs. (|53H56[) phcnomcnologi- 
call y 28 ' 42 . While in the case of a resonant coherent pump, 
one can imagine that the exciton/exciton-polariton dis- 
tribution inherits coherence directly, it is less obvious how 
an incoherent pump can be modelled. A phcnomenologi- 
cal model introduced by Wouters and Carusotto^ 9 - of in- 
coherent pumping has allowed the modelling of the first 
order coherent fraction observed in many experimental 
configurations based on condensatio n 26 ' 60 " —. While the 
GP model cannot model the phase transition during for- 
mation of a condensate and/or supcrfluid, it can offer a 
suitable description of spin currents once spatial coher- 
ence has formed. 

Here we will consider a localized source of the four- 
component indirect exciton system, as corresponds to the 
localized bright spot sources generating exciton conden- 
sates^. We focus our attention on the possible spin po- 
larization textures of coherent excitons propagating away 
from the source. We do not attempt to describe the par- 
tially coherent state within the source, noting that in 
exciton-polariton systems spin currents have been gen- 
erated from both coherent^ 4 - and incoherent^ tightly fo- 
cused spots utilizing the optical spin Hall effect in a 
similar way. We take the polarization state within the 
narrow source area as fixed to a linearly polarized dark 
exciton distribution, assuming that dark excitons have 
lower energy than bright excitons [due to the exchange 
splitting of Eq. (|19p]. The fixed values of the exciton 
wavefunction along the edges of the source area, which 
is assumed circular, acts as a boundary condition and 
effective source for the exciton wavefunction outside the 
"source area. 
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The indirect excitons are known to have very long life- 
time, typically in the range of 10ns to 10/xs^. To model 
the spin currents we thus employ an absorbing boundary 
condition to allow the solution of Eqs. (|53ti56|) in a finite 
area. This can also be seen as allowing a balance between 
source and loss to achieve a steady state of the system (for 
a continuous pump) , where the loss of excitons is caused 
by their travelling outside the area of interest rather than 
their recombination. This is valid for areas of interest up 
to a few (100/im 2 ), bearing in mind that indirect exci- 
tons can easily cover this distance with minimal loss^. 
The steady state solution of the system is independent of 
the initial condition (for our calculations we take = 
in the region outside the source boundary). 

Exciton intensity and polarization distributions calcu- 
lated within the GP approach are shown in Figs. [S] and 
|6l Note that all the plotted quantities are spatially av- 
eraged over 1.5/^m to account for the typical resolution 
of experimental setups^. In principle exciton condensates 
display features on the scale of the de Broglie wavelength, 
A = 27r/|k ex |. Such features are on the sub- micron scale 
and are far beyond experimental resolution. We note also 
that while in the DM approach we could consider exci- 
tons having a fixed radial velocity, in the GP approach 
we necessarily cover the whole range of wavevectors and 
propagation velocities. Also, exciton-exciton interactions 
modify the exciton dispersion. For all these reasons, the 
DM and GP approaches cannot give identical results, in 
principles. The DM approach is suitable for the descrip- 
tion of both coherent and partly coherent exciton gases, 
while, while the GP approach better catches the disper- 
sive propagation features and accounts for nonlinear ef- 
fects. Both approaches are complementary, and it is in- 
structive to compare the results, obtained within these 
two models. 

The exciton density decreases as excitons travel away 
from the source (Fig. Ot). This is not due to exciton 
recombination, which is expected to be very slow, but 
more simply due to the spreading out of excitons in all 
directions. The intensity spread need not be perfectly cir- 
cularly symmetric due to the presence of the spin-orbit 
(Dreselhaus terms) , which can introduce a directional de- 
pendence of the exciton velocity. Fig. \5]p shows the exci- 
ton brightness degree, defined as: 

|vl/ +1 | 2 + |vl/_ 1 | 2 -|vl/ +2 | 2 -|vl/_ 2 | 2 

Pb = 9 9 9 9- (oO) 

l^ + ^+l^^ + l^ + l^l 2 

This quantity represents the degree to which the bright 
exciton density exceeds the dark exciton density. There 
is a conversion of dark to bright excitons as they spread 
out from the source, which can be expected from the 
presence of Dresselhaus coupling terms. 

The polarization distribution is shown in Fig.[6]and can 
be significantly influenced by non-linearity in the system. 
The left-hand plots show the results for negligible non- 
linearity (a = 0; W = 0), which is equivalent to a weak 
pump intensity. Here the polarization distributions are 
qualitatively similar to those calculated in the density 
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FIG. 5: (color online) Spatial distribution of the total exciton 
density (a) and exciton brightness degree (b) [see Eq. (|60|) ] 
in the absence of a magnetic field. The parameters were the 
same as those used for the density matrix DM calculations, 
given in Table [T] with: B = 0; W = 0.2a; m ex was taken 
as 0.21 of the free electron mass. The source area was taken 
circular with a radius of 1/xxn. The images are presented with 
spatial averaging over 1.5/im. The scattering parameter a 
and intensity at the source center were chosen such that the 
interaction energy, a (|*+i| 2 + |*-i| 2 + |*+2| 2 + |*-2| 2 ) = 
1/^eV (being comparable to the other energy scales in the 
system, we are in a nonlinear regime). 



matrix formalism. In analogy to the (intrinsic) spin Hall 
effect^ and the optical spin Hall effec t 10 ' 12 the presence 
of the spin-orbit coupling terms introduces a directional 
dependence of the polarization. The patterns of the po- 
larization degrees divide into quadrants. Some quantita- 
tive differences with the density matrix DM calculations 
appear due to the presence of different wavevectors. 

The right-hand plots show the case of a moderate non- 
linearity, with interaction strength comparable to the 
other energy scales of the system. The most drastic 
effect is on the circular polarization degree, which be- 
comes higher and each quadrant of circular polarization 
divides further giving an eight-lobed pattern to the po- 
larization degree. The interaction terms that we have 
introduced are all spin conserving and it can be noted 
that we have considered the spin isotropic case. Even 
the W nonlinear interaction term, which allows the inter- 
conversion of bright and dark exciton pairs does not ap- 
pear to directly change the spin polarization, conserving 
both circular and linear polarizations upon scattering. 
Still, the nonlinear interaction terms can have a drastic 
effect on the polarization structure. This is because they 
are able to shift (renormalize) the dispersion branches 
in the system. Given that the potential energy of exci- 
tons is fixed by their interaction energy at the source and 
that this energy is converted into kinetic energy at dis- 
tances away from the source, any shifts in the dispersion 
branches can change the wavevector of propagating exci- 
tons. Even if the nonlinear induced shifts of the disper- 
sion branches were not polarization dependent, a change 
in the wavevector of an exciton can allow it to experi- 
ence a different effect from the k-dependent spin-orbit 
coupling terms. In this way, richer structures can appear 
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in the nonlinear regime. 
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FIG. 6: (color online) Spatial distribution of the polarization 
state of excitons in the absence of a magnetic field: horizontal- 
linear polarization degree (a, b), diagonal polarization degree 
(c, d) and circular polarization degree (e, f) for bright ex- 
citons, which corresponds to the near field emission pattern 
of light. The left-hand plots show results in the absence of 
nonlinear interactions (a = 0; W — 0), while the right-hand 
plots show the case of a moderate nonlinearity. The parame- 
ters were the same as in Fig. [5] 
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(61) 



with j = 1,2,3,4 for a = +1,-1, +2,— 2, respectively. 
Experimentally, one can measure the magnetization cur- 
rent associated with the exciton density current. The 
magnetization carried by propagating excitons can be 
found as: 



M, 



PB 

' 2h 



[(9h - 9e) (Pll - P22) + (9h + 9e) (P33 - Pii)} 

(62) 

This expression is obtained having in mind that an elec- 
tron with a spin projection on the z-axis of ±1/2 con- 
tributes to the magnetization projection on the z-axis 
T^5e, and a heavy hole with the spin projection of 
±3/2 contributes to the magnetization T^f-g/i. Hence, 
the magnetization (spin) current produced by the exci- 
tons having a wave- vector k ex will be given by: 



\{9h ~ ge) (/Oll(k eK ) - P22(k e x)) 



2m e „ 

+ (gh + 9e) (P33(k-ex) - paQtex))] (63) 



The total magnetization current in the exciton gas can 
be obtained by integration over all wave-vectors: 



\tot 
JM 



jj U (k ex )dk e 



(64) 



Here A is the area of the sample. This current may be 
detected, for example, by spatially resolved Kerr rotation 
spectroscopy. 




V. EXCITON SPIN CURRENTS 

Consider an exciton state characterized by a wavevec- 
tor k ex and described by the density matrix DM p. Let 
us recall that the elements of this matrix pn, P22, P33, 
P44 are the densities of +1,-1, +2 and -2 spin polarized 
excitons, respectively. The current of each of these den- 
sities is given by a product of the exciton speed and the 



FIG. 7: (color online) Spatial distribution of the exciton spin 
current calculated in arbitrary units using the parameters 
summarized in Table U at k ex — 15.3/im -1 at B — (a) 
and B = 5T (b)-(d) along the z, x and y-axis, respectively. 
The source area was taken circular with a radius of 4/^m. 

Figure [7] shows the spin current density Jm/ (27rr) cal- 
culated for the system of indirect excitons which we con- 
sidered above in the absence of external magnetic field 
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(a) and in the presence of a magnetic field of 5T oriented 
normally to the plane of the structure (b) and along x- 
and y-axes (c,d). In all cases the current intensity de- 
creases as one moves away from the excitation spot, as 
the exciton density decreases inversely with the radius 
r. One can see that the total spin (magnetization) of 
propagating cxcitons experiences oscillations and has a 
strong angular dependence. The spin currents are sup- 
pressed by an in-plane magnetic field, which is not sur- 
prising: jm describes propagation of the normal-to plane 
spin component, which is strongly reduced by an in-plane 
magnetic field. The images show the total spin carried 
by both bright and dark excitons. They do not directly 
correspond to the polarized photoluminescence map for 
two reasons: first, dark excitons do not contribute to 
the photoluminescence; second, the polarization degree 
of photoluminescence does not experience the 1 jr decay 
characteristic of the total spin density. On the other 
hand, the images presented in Figure [7] do correspond to 
the signal of spatially resolved Kerr or Faraday rotation, 
which is sensitive to the normal-to-planc magnetization. 



VI. SPIN CURRENTS IN EXCITON 
CONDENSATES 

The approach formulated above can be extended to 
the description of spin currents in coherent exciton (or 
cxciton-polariton) condensates accounting for particle- 
particle interactions. In this case we need to replace 
the momentum Kk ex by a momentum operator p = 
—ih\J and the diagonal components of the density ma- 
trix DM pn, p 22 , P33, Paa by the exciton densities 
|i/> + i| 2 , | -0 1 1 2 , IVM-2I 2 , |V'-2| 2 , respectively, in the expres- 
sions [HI] and IHll In this case the density currents become: 



a) 



m ex 



(65) 



and the total magnetization current can be expressed as: 
3m = 7^ [(9h - 9e) (niV* + i - xV*-l) 



2m, 



+ (g h + g e ) (*; 2 V* +2 - ** 2 Vtt- a )] 

(66) 



+ (9h+9e) (3+2- j-2)} (67) 

The distribution of the spin density current, j +1 , is 
shown in Fig. |5K. The current density propagates out- 
ward from the source in all directions, decreasing in in- 
tensity. The apparent rotation of the current density is a 
nonlinear effect coming from the interactions in the sys- 
tem. The other spin density currents, j_i, j+ 2 and j+ 2 , 
display a similar behaviour. 

The magnetization current is shown in Fig. [8ji. The 
current is stronger closer to the source, where the intensi- 
ties are stronger. The magnetization current is predicted 
to rotate around the source. 




b) 
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FIG. 8: (color online) Spatial structure of the spin density 
current j+i (a) and the total magnetization current j^-' (b). 
The arrows show the directional dependence of the vector 
fields in space, while the colour code illustrates the intensity. 
The parameters were the same as in Fig. [S] Arbitrary units 
are used for both the spin density and total magnetization 
current. 



One can also introduce the spin conductivity tensor 
linking the components of the density current (|65[) with 
the gradient of potential acting upon each of the exciton 
spin components: 



Ja,l 



l.ni-c-ijj 



(68) 



where I = x,y and m = x,y indicate the in-plane pro- 
jections of the current and potential gradient, respec- 
tively. One can see that c^™ is a 64-component ten- 
sor in the general 2D case. The origin of the poten- 
tial gradient V/7^ jm needs to be discussed separately 
Vf/3,m can originate from the gradient of the quan- 
tum well width, gradient of the barrier height, or it can 
be induced by excitons themselves due to e.g. dipolc- 
dipole repulsion. Indirect excitons have built-in dipole 
moments, the laterally modulated external electric field 
in the z-direction can create an in-plane potential land- 
scape and, in turn,V£/, for them. This was used in 
studies of transport in various electrostatic potential 
landscapes including potential energy gradients^— , cir- 
cuit devices^— , trapsZ^ lattice o 73 i 74 , moving lattices- 
conveyers^, and narrow channel o 71 i 76 ' 77 . 



VII. POLARIZATION CURRENTS 



Spatially resolved measurements of the polarization de- 
grees p c , pi and pd of light emitted by excitons give ac- 
cess to the exciton polarization currents. In terms of 
the density matrix DM formalism, they can be defined 
as products of the exciton speed and the corresponding 
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polarization degree: 





fit 




fit 


Pll 


- P22 




pc = 




Pii 


f P22 


j; (k ea! ) = 


^k ea . 
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Pd = 


hl& ex 


Pll 
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m ex 


Pll 


-f- P22 



(69) 
(70) 
(71) 



The total polarization currents can be obtained inte- 
grating the expressions (|69I71[) over reciprocal space: 



hot 



A 



(27T)' 



/dk eJc , M(kex ) 



(72) 



The polarization currents in an exciton condensate can 



be found from the GP equations (|53M56j) as: 



J/ 



id 



ih 


(*; 1 v* +1 


- q>* 


iV*-0 


Tfl ex 


l*+il 2 


+ |*. 


-1I 2 


ih 
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iV*+i) 
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-ll 2 
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- \&* 
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-ll 2 



(73) 
(74) 
(75) 



The distribution of the polarization currents j c and j; 
are shown in Fig. [5J A striking nonuniform structure 
appears due to the presence of the spin-orbit coupling 
terms. Close to the source spot, there is a strong cir- 
cularly polarized current that rotates around the source. 
This can be attributed to the rotating circular polariza- 
tion degree already observed in Fig. |6] Away from the 
source, the circular polarization current decays, which 
can be expected due to the decay of the spin density cur- 
rent observed in Fig. [5] Along the vertical axis (x = 0), 
a strong circular polarization current remains due to the 
particularly fast change of the circular polarization de- 
gree in this region. The linearly polarized current can 
be stronger further away from the source than at closer 
distances. This is attributed to an increasing linear po- 
larization degree further from the source. In addition, 
one can recall that while spin density currents are gener- 
ally weaker further from the source, there is some com- 
pensation due to the conversion between dark and bright 
excitons (as shown in Fig. [SJd, the bright exciton fraction 
increases further from the source). 



VIII. CONCLUSIONS 

Bosonic spin transport is a young and promising area of 
solid-state physics. The theories of mesoscopic transport 
of charge carriers and quantum transport are among the 
most interesting chapters of modern physics. Substitu- 
tion of fermions by bosons and of a scalar electric charge 
by a spin vector cannot be formally done in these the- 
ories. Basically, all mesoscopic and quantum transport 
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FIG. 9: (color online) Spatial structure of the polarization 
currents j c (a) and j; (b). The arrows show the directional 
dependence of the vector fields in space, while the colour code 
illustrates the intensity. The parameters were the same as in 
Fig- El Arbitrary units are used for the polarization currents. 



effects need to be reconsidered if we speak about electri- 
cally neutral bosonic spin carriers like excitons or exciton- 
polaritons. This is why the area of "spin-optronics" es- 
sentially remains terra incognita. Experimentally, direct 
measurements of indirect excitons and exciton-polaritons 
transport in time-resolved imaging experiments have be- 
come possible in recent years. In this work, we have 
demonstrated that exciton polarization currents are in- 
separably connected with electron and hole spin currents. 
The intensity and direction of exciton polarization cur- 
rents and electron and hole spin currents is governed by 
an interplay of spin-orbit effects, Zecman effects and ex- 
citon exchange effects. In the non-linear regime, the 
pattern of spin currents may also be affected by spin- 
dependent exciton-exciton interactions. 

We have developed two complementary approaches to 
the description of exciton spin currents and textures. The 
DM formalism allows for description of the spin transport 
effects in both classical exciton gases and condensates 
of non-interacting excitons, while the GP equations de- 
scribe propagation of exciton condensates. We predict 
non-trivial topologies of interacting exciton spin in con- 
densates, and suggest tools of their control, such as ex- 
ternal magnetic and electric fields, and source intensity. 
We have demonstrated, that ballistic propagation of ex- 
citons may result in a build up of polarisation patterns, 
which may be observed in near-field photolumincscence 
spectra. 
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